Dynamic behavior of driven interfaces in models with two absorbing states 
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We study the dynamics of an interface (active domain) between different absorbing regions in 
models with two absorbing states in one dimension; probabihstic ceUular automata models and 
interacting monomer-dimer models. These models exhibit a continuous transition from an active 
phase into an absorbing phase, which belongs to the directed Ising (DI) universality class. In the 
active phase, the interface spreads ballistically into the absorbing regions and the interface width 
diverges linearly in time. Approaching the critical point, the spreading velocity of the interface 
vanishes algebraically with a DI critical exponent. Introducing a symmetry-breaking field h that 
prefers one absorbing state over the other drives the interface to move asymmetrically toward the 
unpreferred absorbing region. In Monte Carlo simulations, we find that the spreading velocity of this 
driven interface shows a discontinuous jump at criticality. We explain that this unusual behavior is 
due to a finite relaxation time in the absorbing phase. The crossover behavior from the symmetric 
case (DI class) to the asymmetric case (directed percolation class) is also studied. We find the 
scaling dimension of the symmetry-breaking field yn ~ 1.21(5). 
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Universality classes of models exhibiting a continuous 
phase transition from an active phase into an absorb- 
ing phase with multiple absorbing states are determined 
by the symmetry between the absorbing states If 
the absorbing phase consists of two equivalent absorbing 
states, the phase transition belongs to the directed Ising 
(DI) universality class When the symmetry be- 

tween the two absorbing states is broken, one absorbing 
state becomes completely obsolete and the system crosses 
over to the directed percolation (DP) universality class 
This crossover has been observed in many mod- 
els with two absorbing states, including the probabilistic 
cellular automata (PCA) model |^ and the interacting 
monomer-dimer (IMD) model [Q| in one dimension. 

Dynamic properties of an active domain in the sea of 
the absorbing region are well known for models with a 
single absorbing state, which belong to the DP univer- 
sality class The size of the active domain R in one 
dimension is defined as the distance between the two far- 
thest active sites averaged over surviving samples and 
satisfies the dynamic scaling relation 



At criticality (A = 0), the active domain spreads alge- 
braically in time, R{t, 0) ~ t^^^ with the scaling function 
/(O) being a constant. In the absorbing phase (A < 0), 
the active domain diffuses like a random walker but its 
size converges to a finite value of order of the correlation 
length. So f{x) must scale as {—x)~'^^ in the x —oo 
limit. Therefore the size of the active domain in the long 
time limit scales as 



i?(oo,A)^(-A)- 



(3) 



R{t, A) = ii/^/(At^/'^ii) 



i/^ii 



(1) 



where t is the time, A the reduced external parameter, 
z = the dynamic exponent [ pH , and {i'±) the 

relaxation time (correlation length) exponent. 

In the active phase (A > 0), the active domain spreads 
ballistically with a finite velocity and its size diverges 
linearly in time. So f{x) must scale as x'^^^~^^ in the 
X +00 limit. Therefore the asymptotic value of the 
spreading velocity, v — limt^oo R/t, scales as 



v{A) ~ A''" 



(2) 



and the spreading velocity v = for A < 0. 

For models with two absorbing states, one can consider 
the dynamics of an active domain in two different envi- 
ronments, i.e., in the sea of one absorbing region (defect 
dynamics) and between two different absorbing regions 
(interface dynamics) The defect dynamics describe 

the spreading of a defect (active domain) in a nearly ab- 
sorbing space, while the interface dynamics describe the 
spreading of an interface (active domain) between two 
different absorbing regions. Note that the interface can 
not disappear by itself in contrast to the defect. 

For the symmetric case, the system belongs to the DI 
universality class and the active domain (defect or inter- 
face) behaves similarly in both dynamics. In the active 
phase, the spreading velocity scales as in Eq. (^ with the 
DI critical exponents. At criticality, the size of the ac- 
tive domain grows as R{t, 0) ~ t^^^ with the DI dynamic 
exponent z. In the absorbing phase, the active domain 
splits into two branches, which diffuse like two random 
walkers before they meet and annihilate each other. The 
absorbing domain of the different type from the absorb- 
ing sea emerges between the two active branches. The 
width of each branch should scale as in Eq. (||) with the 
DI exponent but the size of the active domain defined as 
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the distance between the two farthest active sites must 
be of order of the distance between two random walkers. 
So it diverges with the random walk dynamic exponent, 
R ~ tV^Rw with ZRw = 2. 

With a symmetry-breaking field that prefers one ab- 
sorbing state over the other, the system crosses over to 
the DP universality class. The active domain (defect) in 
the sea of the preferred absorbing region behaves in the 
exactly same way as in the defect dynamics of models 
with a single absorbing state, because the probability of 
creating the unpreferred absorbing domain is exponen- 
tially small. The spreading velocity in the active phase 
scales as in Eq. (|2|) with the DP exponents and the de- 
fect size in the absorbing phase as in Eq. (|) with the DP 
exponent. 

However, the dynamics of the active domain (interface) 
between the preferred and the unpreferred absorbing re- 
gion are completely different. The interface is now driven 
into the unpreferred absorbing region by the symmetry- 
breaking field. In the active phase, the interface spreads 
ballistically in both directions. Of course, the interface 
front near the preferred absorbing region (P-interface 
front) moves slower than the interface front near the un- 
preferred absorbing region (JJ-interface front). The un- 
preferred absorbing state is not responsible for the ab- 
sorbing phase transition and its region always tends to 
shrink against the preferred absorbing region. So the dy- 
namic behavior of the U -interface front is always ballistic 
even in the absorbing phase and at criticality. The U- 
interface front velocity, vu, varies smoothly with the ex- 
ternal parameter A and there is no singularity at the crit- 
ical point. The P-interface front behaves like a boundary 
of the active domain in the sea of the preferred absorbing 
region (defect dynamics). The P-interface front velocity, 
vp, scales as in Eq. (^) with the DP exponents in the 
active phase. At criticality, the P-interface front behaves 
in a critical fashion of the DP type. The average spread- 
ing distance of the P-interface front from its initial po- 
sition scales algebraically as t^/^ with the DP dynamic 
exponent z |^]. However, the size of the active domain 
(interface width) is the average distance between the U- 
interface front and the P-interface front, so diverges still 
baUistically even at criticality, R{t, 0) ~ vut. 

As we cross the critical point into the absorbing phase, 
an interesting thing happens. In the absorbing phase, the 
interface still grows ballistically into the unpreferred ab- 
sorbing region with a finite value of vu, but the preferred 
absorbing state dominates the system and tries to con- 
fine the interface within a finite region like the active do- 
main in the DP systems. We expect that the P-interface 
front moves asymptotically with the same velocity and 
the same direction as the [/-interface front in the entire 
absorbing phase. Then vp must have a finite jump at 
the absorbing phase transition and the interface width 
be finite in the long time limit. 

In this paper, we study numerically the dynamics of 
these driven interfaces in the PCA model and the IMD 
model with a symmetry-breaking field. We perform dy- 



namic Monte Carlo simulations to measure the interface 
width and the spreading velocities of the two interface 
fronts. Indeed, we find a discontinuous jump of up at 
criticality and an interesting scaling behavior of the in- 
terface width near the transition. 

The PCA model studied here was originally introduced 
by Grassberger, Krause, and von der Twer ||^. The 
model evolves with elementary rule number 94 in the 
notation of Wolfram except 110 and Oil configura- 
tions, where the central spin 1 flips to with probabil- 
ity p and remains unflipped with probability 1 — p. This 
model has two equivalent absorbing states, i.e., (1010 • • •) 
and (0101 • • •), and exhibits an absorbing phase transi- 
tion that belongs to the DI universality class. We intro- 
duce a symmetry-breaking field h that prefers (1010 • • •) 
over (0101 • ■ •) With probability h, we reject the flip- 
ping attempt of the central spin in the 111 configura- 
tion when it is at an odd-numbered site. As the system 
must go through a 111 configuration right before enter- 
ing into an absorbing state, i.e., (■ ■ ■ 010 111 010 ■ ■ ■) 
(•■• 010101010 •• •), the absorbing state with I's at the 
odd-numbered sites is probabilistically preferable to the 
other for nonzero h. 

For the interface dynamics, we start with a single kink 
(domain wall) between the two different absorbing re- 
gions, i.e., (■•■ 10100101 •■ •), where the sharp interface 
(active domain) with zero width is placed in the middle 
of 00. Then the system is updated in parallel following 
the automata rule. In contrast to the defect dynamics, 
the system never enters an absorbing state. We measure 
the positions of the P-interface front and the U -interface 
front, averaged over 2 x 10'^ samples up to typically 10^ 
time steps for various values of p and h. 

Fig. 1 shows typical evolutions of the interface for (a) 
A > 0, (b) A = 0, and (c) A < where A = p - and 
the critical probability Pc = 0.3908(5) for h = 0.5 [^. 
Evolutions of the average positions of the P- and U- 
interface fronts are shown in (d). The [/-interface front 
always moves ballistically and its velocity vu changes 
smoothly with A. However, the P-interface front shows 
an abrupt change in its dynamics at criticality. It moves 
ballistically in the active phase and its velocity vp van- 
ishes algebraically at criticality. The P-interface front 
still moves toward the preferred absorbing region at crit- 
icality and its average distance from the initial position 
diverges as t^^^ with the DP dynamic exponent z. In the 
absorbing phase, the P-interface front turns around and 
moves ballistically in the opposite direction. In the long 
time limit, the P-intcrface front moves with the same ve- 
locity and the same direction as the [/-interface front for 
A < 0. 

In Fig. 2, we plot vp and vu versus A for h — and 
0.5. The critical probability Pc = 0.125(2) for the sym- 
metric case {h = 0) ||l^. Indeed, we find that there is 
a discontinuous jump of up at criticality for the asym- 
metric case {h = 0.5) and both velocities coincide in the 
absorbing phase. In the active phase, vp vanishes con- 
tinuously approaching the critical point like in the sym- 
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metric case. Log-log plots of vp versus A near criticality 
in the active phase are shown in Fig. 3. From the slopes, 
we estimate v\\ — v±^ as 1.35(10) for h ^ and 0.65(3) 
for h — 0.5, which agree well with the DI (1.39) and DP 
(0.637) value, respectively. We also run dynamic Monte 
Carlo simulations for the IMD model and find the similar 
results, i.e., 1.25(15) for the symmetric case and 0.62(3) 
for the asymmetric case. 

In the absorbing phase, the average distance between 
the P-interface front and the t/-interface front (inter- 
face width) is found to saturate to a finite value, Rs — 
R{oo,A). However, Rs does not scale as in Eq.(i be- 
cause this length scale is not proportional to the corre- 
lation length of the system. As one can see in Fig. 1(c), 
this length does not measure the typical size of the DP 
active domain but the distance from the unpreferred ab- 
sorbing region to the farthest active site in the preferred 
absorbing region. As the ?7-interface front always moves 
ballistically into the unpreferred absorbing region, this 
distance should be proportional to the time scale in the 
system which measures the average lifetime of the tree- 
like active DP clusters grown out of the [/-interface front. 
However, small clusters do not contribute effectively to 
the average lifetime due to the shading effect by big clus- 
ters. This effect complicates the scaling behavior of Rs- 

A simple-minded scaling theory for R{t, A) in the 
asymmetric interface dynamics may be given as 

R{t,A)=tg{At'^''«), (4) 

where is the relaxation time exponent for the DP uni- 
versality class. At A = 0, g{0) is a constant, so R{t, 0) 
diverges linearly in time, which is correct for the asym- 
metric case. For A < 0, R{t, A) should saturate in the 
long time limit, so g{x) must scale as {—x)~'^^^ in the 
X — > — oo limit. Therefore we expect that Rg scales as 
(-A)-''" with the DP exponent v\\ =i 1.733. The log- 
log plot of Rs versus |A| shows a fairly good straight 
line, which seems to support a simple power-law scaling 
of Rs (see Fig. 4). However, our estimate for the scal- 
ing exponent from its slope is well over the above DP 
value, i.e., 2.00(5) for the PCA model and 1.95(10) for 
the IMD model. We find that this discrepency is due to 
the shading effect. A careful analysis incorporating the 
shading effect suggests that there should be a logarithmic 
correction in the scaling theory, i.e., Rs ^ lAj^^ii ln|A| 
0. The log-log plot of Rs/\\n\A\\ versus |A| (Fig. 4) 
also shows a fairly good straight line. From the slope, 
we estimate the scaling exponent v^^ = 1.75(5) for the 
PCA model and 1.65(10) for the IMD model, which arc 
in good agreement with the DP value. 

Finally, we study the crossover behavior from the DI 
class to the DP class. The operator associated with the 
symmetry-breaking field h must be relevant at the DI 
critical point, so the scaling dimension of this crossover 
operator must be positive and may be an indepen- 
dent DI critical exponent. We obtain yh numerically by 
measuring the interface front velocities at small values 



of h along the p = line, where p^ = 0.125 is the DI 
critical probability at /i = 0. At finite values of h along 
the Pc — Pc line, the system becomes absorbing and the 
two interface front velocities, vu and vp, coincide and 
become finite. 

Consider the crossover scaling relation near the DI crit- 
ical point for the average position of the unfavorred in- 
terface front Xij measured from its initial position as 
Xu{t,h) = bXu{b-''°'t,by'-h) where zdi is the DI dy- 
namic exponent (1.75) and b an arbitrary scaling factor. 
With b = we have Xu = ti/^"^$(W?^''/^°0- For 

the symmetric case {h — 0), Xu scales as t^l^i^' with 
$(0) being a constant. For ft, > 0, Xu increases linearly 
in time (finite Vu), so <^{x) ~ x'^"^^^)/^'' in the a; ^ oo 
limit. Therefore, vjj ^ h'^ where n — {zjji — l)/yh. 

Log-log plots of Vjj versus h at p — p*^ are shown in 
Fig. 5. From the slope, we estimate k — 0.62(2) and 
hence obtain yh — 1.21(5). Interestingly, this value is 
very close to the value of the scaling dimension of the 
operator associated with a roughening degree of freedom 
at the DI critical point ||l^ . The conventional definition 
of the crossover exponent cj) is given as the ratio of the 
two scaling dimensions, i.e., ip = yu/yA where j/a is the 
DI scaling dimension of the temperature-like operator A 
and j/A = Using this relation, we find (f) — 2.24(10) 

which is consistent with the recent result for general- 
ized monomer-monomer models studied by Bassler and 
Browne Q. 

In summary, we study numerically the dynamics of 
the driven interfaces in the PCA and IMD model that 
have two asymmetric absorbing states. We find that the 
spreading velocity of the driven interface shows a dis- 
continuous jump at criticality due to the finite interface 
width in the absorbing phase. The interface width di- 
verges in a nontrivial manner, approaching the critical- 
ity. We find that our numerical data are consistent with 
a recent scaling theory taking into account the shading 
effect of big active clusters over small ones [|4| . 
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FIG. 1. Typical evolutions of the asymmetric interface 
dynamics for (a) A = 0.05, (b) A = 0, and (c) A = -0.05 
with h — 0.5. The preferred (unpreferred) absorbing region 
is shown in black (grey) and the active sites are represented 
by white pixels. Evolutions of the interface fronts averaged 
over many samples are shown in (d): dotted lines for the ac- 
tive phase, solid lines at criticality, and dashed lines for the 
absorbing phase. 
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FIG. 2. Interface front velocities versus A. Filled circles 
and open boxes represent the interface front velocities for the 

symmetric case. Filled up-triangles represent the P-interface 
front velocity vp and open down-triangles the {/-interface 
front velocity vu for the asymmetric case {h = 0.5). Lines 
between data points are guides to the eye. 




FIG. 3. Log-log plots of Vp versus A. Filled circles are 

for the symmetric case and filled triangles for the asymmetric 
case {h = 0.5). The solid lines are of slope 1.35 and 0.65, 
respectively. 
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